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ABSTRACT 

We present a spectral analysis of the narrow-line Seyfert 1 galaxy I Zwicky 1 , focusing 
on the characteristics of the ionised absorbers as observed with XMM-Newton in 2005. 
The soft X-ray spectrum shows absorption by two components of ionised gas with a 
similar column density (Nh ~ 10 21 cm -2 ) and ionisation parameters log£ ~ and 2.5. 
Comparing this observation with a 2002 XMM-Newton data set, we see a clear anti- 
correlation between the X-ray ionisation parameter £x and the 0.1 — 10 keV luminosity. 
Viable explanations for this effect include transient clouds or filaments crossing the 
line of sight in a complex geometry or a gas observed in non-equilibrium. The outflow 
velocity of the X-ray low-ionisation absorber is consistent with the outflow of the UV 
absorber detected in a past Hubble Space Telescope observation. In addition, the ionic 
column densities of C iv and N V derived from the X-ray model are consistent with 
the UV values. This suggests that the low-ionisation outflowing gas may survive for 
many years, despite large changes in flux, and that there is a tight connection between 
the X-ray and UV absorbers that can only be confirmed with a simultaneous UV and 
X-ray observation. 

Key words: Galaxies: individual: IZwl - Galaxies: Seyfert - quasars: absorption 
lines Ultraviolet: galaxies X-rays: galaxies 



1 INTRODUCTION 

About 50-70% of active galactic nuclei (AGN) display ab- 
sorption by ph o toionised gas in th eir X-ray spectra (e.g. 
iRevnoldsl Il997l ; iGeorge et ail Il998j ). With the advent of 
high energy resolution instruments, the study of the bright- 
est sources has a llowed us to detail the characteristics 
of such gas (e.g. iKaastra et ail |2002l: iKaspi et al.l 120021 : 
IKraemer et al.ll2005l ; iTurner et al.ll2005h . Multiple compo- 
nents are often detected, differing in column density, ionisa- 
tion and outflow velocity. At least some of the lower ion- 
isation X-ray components most proba bly share the same 
location as the UV c omponents (e.g. ITurner et al.l 120051: 
ICostantini et al] l2007f ) as first suggested by Mathur et al.l 
l| 19951 ). Variability studies point out the complex response 
of the ionised absorbers (also called warm absorbers) to 
central-source flux variations. A linear response of the 



ionisation parameteiQ of the gas component to the con- 
tinuum change pr ovides a constraint on the distance to 
the absorber (e.g . iNetzer et all 120021 : ITurner et ail 120051 : 
lOtani et al.lll99a ). A linear response gives a relatively in- 
tuitive picture of the geometry of the absorber as a thin 
slab of material in thermal equilibrium (either in clouds or 
filaments) exposed to th e ionising flux. In ot her variability 
studies (e.g. NGC 4151. IKraemer et 81.1120051 ) the gas does 
not obey this simple picture, suggesting instead that gas 
outflows often change direction along the line of sight in a 
unpredictable way. 

I Zwicky 1 (I Zw 1) i s a well-studied narrow-li ne Seyfert 1 
galaxy (NLS1; see IVestergaard Wilkes! 1200 ll . and refer- 
ences therein) loca ted at redshift z=0.06112±0.000015 
jCondon et all Il985l ). Since the time of ROSAT observa- 
tions, I Zw 1 has been known to host intrinsic absorption 
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1 The ionisation parameter is defined as § = L/nr 2 , where L is 
the 1-1000 Ryd ionising luminosity, n and r are the density and 
the distance to the ionised gas. 
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in the X-ray band jBoller et al.lll996l ; lLawrence et al.|ll997l ). 
From moderate-resolution ASCA data, Leighlvl ( 19991 ') could 
distinguish between a neutral component (Nn ~ 2.6 x 
10 20 cm -2 ) and an ionised absorber, typically highlighted, 
in lower resolution spectra, by the Ovu and Oviii absorp- 
tion edges. A sim ilar spectral complexity was reported by 
iGallo et all (|2004l ), using XMM-Newton data. The presence 
of fast (uout ~ — 1870 kms -1 ) outflowing ionised gas was 
observed in the UV spectrum of I Zw 1 (using HST-FOS; 
lLaor et "al]|l997l ; ICrenshaw et al.lll999l ). through the analy- 
sis of the N v and C IV absorption troughs. FUSE observed 
the O VI absorption troughs, which displayed multiple com- 
ponents i|Krissll2002r i. 

This paper is one in a series devoted to a long XMM- 
Newton observation of I Zw 1. In the present study we focus 
on the soft X-ray absorbers, discussing the physical picture 
returned by the long-term variability of this gas and its con- 
nection with the UV absorbers. Other works deal with the 
continuum emission and the variable iron Ka line (Gallo et 
al. 2007a, hereinafter Paper I), and with the spectral vari- 
ability analysis (Gallo et al. 2007b). 

The paper is organised as follows. In Sect. 2 the main 
features of the data reduction and the EPIC-pn analysis 
(described in detail in Paper I) are summarised. Sect. 3 is 
devoted to the analysis of the RGS data and the modelling 
of the spectral features in IZwl. In Sect. 4 we discuss our 
results and in Sect. 5 we present the conclusions. The cosmo- 
logical parameters used are #o=70 kms^Mpc" 1 , O m =0.3, 
and f2A=0.7. The quoted errors are at 90% confidence for 
one interesting parameter, unless otherwise stated. The as- 



and hard (2-10 keV) energy bands, respectively. The flux 
changes after ~27ks since the start of the observation and 
happens on a dynamical time scale of ~5 ks. The flux varia- 
tions in the soft and hard X-ray bands show similar trends 
(see Fig. 3 of Paper I). We will study the effect of the flux 
change on the RGS spectrum in Sect. [3] 
The data from the EPIC-pn camera were used to constrain 
the broad-band spectrum. The modelling of the continuum is 
described in full detail in Paper I. For the purpose of study- 
ing any narrow absorption/emission features in the RGS 
band, the detailed broad-band continuum shape is not cru- 
cial. For simplicity, here we describe the continuum with a 
broken power-law, which fits the data well (Tab[T}. However 
in Sect. 13. D we test our results for different continua. 
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sumed Galactic c olumn density is 
l|Elvis et al.lll989l ). 



2 THE DATA ANALYSIS 
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3 THE SOFT ENERGY SPECTRUM 

The RGS spectrum shows the signatures of absorption, 
pecially around the O VII and O viii absorption edges 
order to fit a photoionised absorber, we used the SPEX 
model XABS, which interpolates on a fine grid the TVh and 
log ^ values. The ionisation balance is constructed from the 
spectral energy distribution (SED) sp ecific for the source 
(solid curve in Fig. [!}, using Cloudy (|Ferlandll2004 ). Here 
and in the following analysis, the SED definition takes ad- 
vantage of the simultaneous EPIC-pn and OM data. The 
OM data were collected in fast-mode (see Paper 1) in the 
U (3440A), UVW1 (2910A), UVW2 (2120A) and UVM2 
(2310A) filters. We considered the average flux for each fil- 
ter. By interpolation, we obtained the flux at 2500 A that we 
used in constructing the SED. Interpolating the OM points 
was mandatory in order to compare consistently the ionisa- 
tion balance obtained with this SED with the one obtained 
from a previous XMM-Newton data set for which only one 
OM filter was available (Sect.E2J. An adopted SED will al- 
ways have a degree of uncertainty, when dealing with limited 
multiwavelength data. When comparing SEDs of different 
flux states of the same object it is important to consistently 
compare the resulting ionisation balances. The SED of 2002 
was likely similar to 2005 in the UV regime, but in 2002 
only one OM filter was available. This would make any "re- 
construction" of the 2002 SED highly arbitrary. Therefore 
the only way to consistently compare the SEDs at the two 
epochs is to interpolate the 2005 OM points. Beyond the 
optical band, we used a standard AGN continuum, char- 
acterised by relatively low emission at longer wavelengths, 
as defined in Cloudy. The host galaxy of I Zw 1 is indeed 
very luminous both in the far-infrared and in the infrared 
l|Halpern Okelll987l ). In the long-wavelength band it is im- 
portant to avoid dust contamination of the nuclear emission 
by the host galaxy. In general, an overestimation of the in- 
frared radiation may le ad to an incorrect io nisation balance 
determination (see e.g. lFerland et al1l2002t ). 
Tunable parameters of the XABS model, in addition to Nu 
and log£, are the root mean squared width of the absorp- 
tion lines, the outflow velocity, the line-of-sight covering fac- 
tor of the gas and the elemental abundances. We fixed the 

|http: //xmm.vilspa.esa.es/external /xmm_u S er_support /user S group/2006^ lgyWHefe ctor to „ un ty and we assumed solar abundances 
shtml. ijAnders Grevesse 1989) throughout the analysis. 

3 http://www.sron.n1/divisions/hea/spox/version2.0/ In Tab. Q] we list the best-fit parameters for the absorbers 

release/index. html in I Zw 1 as measured by the RGS and, for consistency, with 



1027) 
cam eras MOS 
, the high 

resolution spectrometer R GS (den Herder et alj|200ll ) and 
the optical monitor (OM, iMason et al.ll200ll ) were success- 
fully gathered. After filtering high-background intervals, the 
net exposure time was ~ 58 ks for EPIC-pn and ~ 84 ks for 
RGS. The event files for each instrument were created with 
the newest SAS pipeline (ver. 7.0). This version of SAS ben- 
efits from improved calibrations of RGS time-dependent ef- 
fective area, absolute RGS flux and bad-pixel location^. The 
discrepancies between EPIC-pn and RGS are of the order of 
10%, similar as for EPIC-MOS when compared to EPIC- 
pn. 

The spectral analysis of the soft X-ray spectrum was car- 
ried out using the fitting package SPEX0 (ver. 2.0). We re- 
binned the RGS data to reach a signal-to- noise ratio (S/N) 
of ~9 over the whole energy band. During the observation 
the source shows two different average flux states, which dif- 
fer approximately by 44% and 36% in the soft (0.3-2 keV) 
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the lower resolution EPIC-pn data. The model parameters 
are always consistent between RGS and EPIC-pn. However, 
the higher energy resolution makes RGS more reliable for 
the modelling of narrow absorption features, and we will 
consider the physical parameters constrained with the RGS 
data for the discussion. 

A warm absorber (hereinafter labelled 05a), with a column 
density Ah ~ 1.3 x 10 21 cm -2 and an ionisation parame- 
ter log ^ ~ is highlighted in the RGS data especially by 
the blend of the O VII K absorption edge and the iron un- 
resolved transition array (UTA). The gas is found to be 
outflowing at ~1700kms - . The widths of the lines are 
loosely constrained (a — 36 ± 20 km s -1 ), but still consistent 
with the more accurate FUSE measurement of 30±5kms -1 
(G. Kriss, private communication). Additional ionised ab- 
sorption is required to fit the region below the Oviii edge 
(A = 14.22A). The improvement of the fit after the inclu- 
sion of this second absorber is Ax 2 /A^=27/2, where v is 
the number of degrees of freedom. The second gas com- 
ponent (hereinafter 05b) has a higher ionisation parame- 
ter (log£ ~ 2.6, Tab. [TJ but a column density similar to 
the 05a component. We fixed the widths of the lines to the 
value for the lower ionisation component. The absorption 
lines of the higher ionisation component lie in a spectral re- 
gion of lower resolution, and a possible outflow velocity is 
hard to quantify. The upper limit on the outflow velocity of 
the 05b component is v ou t < 800 km s - . Finally, a neutral 
component, identified by an O I edge at the redshift of the 
source and with A^h ~ 1.5 x 10 20 cm -2 , is required by the 
data. The addition of this component improves the fit by 
Ax 2 /Aj/=19/1 (where the column density is the additional 
parameter). In Fig. [2] we plot the transmission of the neu- 
tral and warm gas components (upper panel), the observed 
RGS spectrum with the best-fit model (middle panel) and 
the residuals from the best fit (lower panel). 
The above analysis refers to the total RGS spectrum. How- 
ever (Sect. [2}, the soft and hard energy bands undergo vari- 
ations during the XMM-Newton pointing. Using the RGS 
data we further verified that only continuum changes are re- 
sponsible for the overall variations in the source flux and not 
opacity variations of the ionised gas. We divided the RGS 
data into two time segments, corresponding to the higher 
(~30ks long) and lower state (~40ks). The higher state 
best fit does not vary significantly with respect to the aver- 
age spectral fit, except for the normalisation and slope of the 
power law. We then used the higher flux state best-fit model 
(see above) as a template for the lower flux state and found 
that the neutral and warm absorbers do not vary within the 
uncertainties. 

3.1 Other continuum models 

In the present paper the continuum is represented by a sim- 
ple broken power law. A discussion of the continuum models 
and their physical implications is given in Paper I. Here we 
test different continuum models only in order to assess the 
robustness of the warm-absorber parameters. We follow the 
same approach as in Sect. [3j the broad-band continuum is 
established using EPIC-pn and then applied to the RGS, in 
order to tune the warm absorbers' parameters. 
First we used a power law plus modified black-body emis- 
sion and neutral absorption (both from the host galaxy and 
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Figure 1. The adopted SEDs for the 05 observation (solid line) 
and 02 observation (dashed line). The X-ray shape of the SED is 
based on the EPIC-pn data, while the point at 2500A is derived 
by interpolation of the OM data points. The 05 OM points (open 
diamonds around 10 15 Hz) were interpolated to the same wave- 
length measurement of 02, in order to compare consistently the 
05 SED with the 02 one, for which data from only one OM filter 
was available. 



Table 1. 2005 best-fit model parameters for RGS and PN. Fluxes 
are in units of 10 — 12 ergs -1 cm -2 , unabsorbed luminosity in 
units of 10 43 ergs -1 . Component 1: broken power law, 2: neu- 
tral absorber, 3,4: warm absorbers. All emission components are 
modified by Galactic photoelectric absorption. 



comp. 


par am. 


R1+R2 


PN 


1 


r 2 

-Ebrcak (keV) 


2.5±0.1 


2.58±0.05 
2.21±0.05 
2.0±0.1 


2 


N H (10 20 cm -2 ) 


1.7±1 


1.5±0.4 


3 


N H (10 20 cm -2 ) 


13±3 


9±1 




log £ 


0.05±0.16 


0.16±0.14 




v out (kms -1 ) 


-1700±400 


-1700 fix. 


4 


Nn (10 20 cm -2 ) 


13±5 


8±3 




log? 


2.6±0.3 


2.0±0.4 




X> 


411/358 


530/513 




F2-10kcV 




4.86±0.06 




Fo.5-2 keV 


5.0±0.3 


5.13±0.07 




L 2 -10 keV 




4.50±0.06 




Lo.5-2 keV 


7.0±0.4 


6.10±0.08 



from the Milky Way). In the black-bo dy model the effect s 
of Compton scattering are considered (|Kaastra Bardll~989r ). 
A black-body temperature of ~0.13keV is required to fit 
the low energies. Without any warm absorber the result- 
ing reduced \ 2 (x 2 l v ) ls ~l-32. The addition of one warm 
absorber (two additional free parameters) improves the fit 
by Ax 2 =93. A second component improves the fit by an- 
other Ax 2 /Af = 23/2. The final model has a higher black 
body temperature, kT ~0.2keV (see Paper I for a discus- 
sion). The total black-body luminosity is however very low 
(~ 2.7 x 10 42 ergs -1 ) compared to the dominant power law 
component (~ 1.3 x 10 43 ergs -1 ) in the 0.3-10 keV band. 
The warm absorbers parameters do not differ within the er- 
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Figure 2. Upper panel: transmission models (convolved with the data resolution) of the neutral absorber (dashed line) and the ionised 
components for the 05 spectrum. Component 05a is displayed with the light line and 05b with the dark line. The most relevant transitions 
are also labelled. Middle panel: RGS data together with the best fit. Lower panel: residuals from the best-fit model. 



rors from the ones listed in Tab. [T] 

Next, we tested a relativistically blurred reflection model 
jRoss Fabianll2005l ). We fixed the parameters to the model 
values determined in Paper I, i.e. the disc inclination to 51°, 
its emissivity a = 2.5 and inner radius to 5 gravitational 
radii. We assumed solar abundances for the disc material. A 
fit without any ionised absorption results in a % 2 l v °f 1-56. 
This value decreases to 1.21, adding one warm absorber. A 
second warm absorber component again improves the fit, as 
in the previous cases (Ax 2 /A^ = 20/2). Note that for this 
exercise we do not attempt any fit of the Fe Kq line. The 
amount of reflected light with respect to the total luminosity 
budget is modest (few percent), resulting in a weak soft ex- 
cess. We conclude that for the purpose of our study a simple 
broken-power law fit is acceptable. 



3.2 The recent history of the absorbed X-ray 
spectrum 

I Zw 1 was observed by XMM-Newto n for the first tim e on 
June 22, 2002 (rev. 0464) for 20 ks (|Gallo et alj|2004 ). At 
the time the source 2-10 keV flux was about 1.7 times higher 
than in the present observation and showed remarkable spec- 
tral complexity at low energies. In Fig. we over-plot the 
two EPIC-pn data sets, but with the 2005 spectrum shifted 
up in order to best match the continuum during 2002. A 
qualitative inspection of the figure shows that the soft ener- 
gies underwent the most dramatic changes between the two 
epochs (we refer to the 2005 and 2002 data sets as 05 and 




1 

Energy (keV) 



Figure 3. EPIC-pn spectra taken in 02 (light line) and 05 (dark 
line). The spectrum of 05 is multiplied by 1.7 to best match the 
hard continuum during 2002 and emphasises the low energy com- 
plex variability. 



02, respectively). We reprocessed the 02 RGS and EPIC-pn 
data using SAS 7.0 and we focused the analysis mostly on 
the warm absorber. The RGS data were rebinned to reach 
a S/N ratio ~ 9 over most of the band. The best fit shows 
that the underlying continuum changed between 02 and 05 
mostly in flux rather than shape, which is here again de- 
scribed by a simple broken power law (Tab.[2)|. The shape 
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Figure 4. Upper panel: transmission models (convolved with the data resolution) of the neutral absorber (dashed line) and the ionised 
components for the 02 spectrum. Component 02a is displayed with the light line and 02b with the dark line. The most relevant transitions 
are also labelled. Middle panel: RGS data together with the best fit. Lower panel: residuals from the best fit. 



of the ionising continuum was again constructed using the 
unabsorbed EPIC-pn continuum and the fluxes recorded in 
optical/UV by OM (dashed line in Fig. [TJ. In 02 only the 
OM-UVW2 filter was used. Similarly to the 05 data, in ad- 
dition to the neutral gas, we find evidence for two ionised 
absorption components. The corresponding ionisation pa- 
rameters are log£ ~ —0.9 and log£ ~ 1.6. We identify 
these two warm absorber components as 02a and 02b, re- 
spectively (Fig. [4} . Interestingly, we find evidence that the 
02a component outflows with a velocity similar to the 05a 
component (« ou t ~ 1800 kms" 1 ). The width of the lines is 
a = 67 ± 50 kms - . Going further back into the recent his- 
tory of the X-ray absorbers o f IZwl, we r eanalysed the 
ASCA data obtained in 1995 (|Leighlvi [l999l ) . The contin- 
uum is once again well fitted by a broken power law with 
parameters comparable to the most recent 05 XMM-Newton 
spectra: T soft = 2.58 ± 0.07, r hard = 2.35 ± 0.07 and a break 
energy at 2.0 ± 0.8 keV. We confirm the presence of a neu- 
tral gas phase with column density Nn = 4±1 X 10 20 cm" 2 , 
consistent with the values found with XMM-Newton. 
For the ASCA data, simultaneous UV observations were un- 
available; therefore a true SED could not be constructed. 
The ionising SED was then built using the ASCA contin- 
uum plus the same a O x0 = —1.25 found in 05. We used this 
SED only to obtain a best fit that consistently accounted for 
the absorption features and then only considered the ionic 



column densities (and not the ionisation parameter) to com- 
pare with our XMM-Newton results. 

As the ionisation parameter is dependent on the ionisation 
balance (i.e. the SED), a good way to compare warm ab- 
sorbers produced by different (or incomplete) SEDs, is to 
evaluate the different elemental ionic column densities. In 
particular, iron is the most sensitive to changes in the ioni- 
sation of the gas (e.g. iBehar et al]|200ll ). In Fig. [SJ we plot 
the ionic column densities as a function of the ionisation 
stage of iron for each of the X-ray ionisation components 
detected in this source. The column densities are obtained 
from the XABS best-fit models. The relative errors are sim- 
ply the scaled errors from the total hydrogen column den- 
sities of the XABS ionisation components. The 02 compo- 
nents are well separated: 02a peaks at Fev, while 02b peaks 
at Fexiv-Fexvi. The 05a component has a relatively lower 
total column density (Tab. [TJ and shares some ions with 
the 02a components. The 05b component displays the high- 
est ionised iron transitions. The ASCA component (labelled 
95) is more in the high ionisation region. 
A possible second component in the ASCA spectrum would 
be difficult to detect. We verified that a higher ionisation 
component, with a column density of the order of 10 21 cm" 2 
would produce mostly lines (see Fig. [2j [4| , which would be 
undetectable in the ASCA spectrum. A much less ionised gas 
peaking, say at Feni-Fev, would be easily confused with the 
neutral absorber. 



4 a x is the spectr al index between 2500 A and 2 keV 
llTananbaum et al.lfl979h . 
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Figure 5. Ionic column density as a function of the ionisation 
level of iron in the warm absorber components detected in recent 
years by XMM- Newton (02a, 02b, 05a and 05b) and ASCA (95). 



Table 2. 2002 best-fit model parameters for RGS and PN. Fluxes 
are in units of 10 — 12 ergs -1 cm -2 , unabsorbed luminosity in 
units of 10 43 ergs -1 . Component 1: broken power law, 2:neu- 
tral absorber, 3,4: warm absorbers. All emission components are 
modified by Galactic photoelectric absorption. 
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Figure 6. Stability curve for the 2002 (solid line) and 2005 
(dashed line) observations. The dark circles represent the two 
ionised components detected in 2002, while the light circles rep- 
resents the 2005 components. None of the components share the 
same pressure ionisation parameter. 
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Figure 7. The log of the X-ray ionisation parameter, §Xj against 
the 0.1-10 keV intrinsic luminosity for the 02 (dark dots) and 05 
(light dots) components. 



4 DISCUSSION 

4.1 The X-ray warm absorber 

The 02 and 05 spectra clearly show absorption by ionised 
gas. There were two detected components in 02, with ioni- 
sation parameters log£= -0.9 and 1.6. The 05 spectrum was 
also well fitted with two warm absorber components, show- 
ing absorption with higher ionisation (log£ ~ and 2.6). In 
Fig. [6] we plot the gas temperature (T) against the ionisa- 
tion pressure parameter S for the 02 and 05 SEDs (drawn 
with the solid and dashed lines, respectively). The ionisa- 
tion pressure parameter is defined as 3 = L/4irr 2 cp, where 
L is the ionising luminosity, r the gas distance from the 
source, and p the gas pressure. The dark dots mark the po- 
sition of the 02 warm absorbers on the stability curve rel- 
ative to the SED at that epoch, while the light dots are 
the components detected in 05. Our curves show a wide 
interval of logH in which the absorbers exist in a cold 



regime. This is where the "cold" components 02a and 05a 
lie. The ascending part (log 3 > 0.5) of th e curve does not 
prese nt the typical turnover S-shape (e.g. ICrenshaw et al.l 
Il999h . T his is because of the relatively steep X -ray SED 
of I Zw 1 (jKrolik et al.lll98ll ; iGuilbert et ai]|l983l ). We note 
that both the two 02 and two 05 components lie far apart 
from each other and are not in pressure equilibrium. This 
beha viour has been observed in other Seyfert 1 g alaxies 
(e.g.. ISteenbrugge et al1l2005l ; ICostantini et all 120071 ). This 
suggests that cond itions other than pressure equilibrium 
l|Krolik KrissHl995l ) should be present to make the two ab- 
sorbers co-exist, such as magnetic confinement , as alread y 
suggested for the broad line region clouds (e.g. |ReesJ|l987|). 
The lower ionisation a components show an outflow velocity 
of roughly 1700-1800 kms" 1 in 02 and 05. Therefore there 
are no strong velocity changes apparent over time, although 
the constraints on velocity changes are limited. 
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4.2 Long-term changes in the warm absorbers 

The X-ray ions are mostly influenced by the 0.1 -10 keV 
radia tion, rather than by the UV photons (e.g. iNetzer] 
1996). Therefore in the discussion below we shall use 
the X-ray ionisation parameter £x> obtained simply by 
scaling do wn the ionising lu minosity to the 0.1-10 keV 
band (e.g. iGeorge et al.l 19981 ). This calculation can be 
done only for the XMM-Newton data sets, for which the 
complete ionising SEDs were available. The values for £x 
then become -1.96 ± 0.2 and 0.43 ± 0.22 for the 02a and 
02b components, respectively. The 05a and 05b values are 
-1.16 ± 0.16 and 1.43 ± 0.3, respectively. In Fig. [7] £x of 
both the 02 and 05 components are plotted against the 
0.1-10 keV luminosity. Note that in general the ionisation 
parameter value is SED dependent. For the 02a and 05a 
components the corresponding stabili ty curves (which are 
mostly dependent on the SED shape, iBottorff et al.l |2000| ) 
are indistinguishable (Fig. [6]). The stability curves start 
to differ slightly at higher values of H a (,/T, therefore 
we compare the 02b and 05b results with this caveat in 
mind. We note a significant anti-correlation between the 

0. 1 — lOkeV luminosity and the X-ray ionisation parameter. 
For a twofold decrease in luminosity between 02 and 05, the 
low ionisation (a components) £x increases by a factor of 
~ 6.2. For the 02b and 05b components, the increase factor 
is ~ 10. We may consider the a and b absorbers of 02 and 05 
to be low- and high-ionisation components. If the observed 
ionisation comes from gas in equilibrium, we should invoke 
either a decrease in density (i.e. the gas thickness should 
increase, as Ah did not change dramatically) between 02 
and 05, or movement of the gas closer to the central source. 
The latter seems an unlikely possibility as blue-shift of the 
a components is observed in 02 and 05. If the gas consists of 
discrete outflowing narrow filaments or clouds, then we do 
not expect to see the same structure along our line of sight 
if the distance Ar travelled by this gas is much larger than 
its characteristic size d. For a rest-frame time difference of 
~1057 days between both observations and the measured 
outflow velocity of 1700 kms -1 , Ar ~ 1.5 x 10 16 cm. 

For the following calculation we consider the case Ar > d, 

1. e. the size of the cloud is comparable to the gas travel 
path. The absolute distance of the centre of the cloud 
r from the source should be much larger than the sizes 
considered here, therefore we also assume that the change 
in r is not a critical term. We then compute 



£05 

£02 



where nd — Ah- The observed £x ratio between 05 and 02 
is 7 for the a components. The X-ray luminosity ratio is 0.5. 
With a simple calculation we obtain dos/do2 = 7.6. Keep- 
ing in mind the assumption Ar > do2 ~ 1.5 x 10 16 cm, 
for the a components we obtain no2 < 1.6 x 10 5 cm -3 and 
nos < 1.1 x 10 4 cm -3 . These are very low values com- 
pared to the estimated l imits for the warm-absorber den- 
sities 2 x 10 6 - 10 9 cm -3 (|Netzer et al.ll2002l ; IKraemer et al.1 
2002), therefore it seems unlikely that we are observing the 
same portion of the gas in both 02 and 05. 
A scenario that could reproduce the observed anti- 
corr elation sees the gas co nstantly in a non-equilibrium state 
(e.g. INicastro" et all 1 19991 ) because, for example, of a delay 
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Figure 8. Column density of hydrogen (Hl+Hn), Nv and Civ 
predicted by the X-ray absorbers at different epochs. The dashed 
vertical lines show the ra nge of TVh inferr ed from the HST mea- 
surements taken in 1994 llLaor et al.lll997l) . 



in the gas response. Such a delay may be due to the gas hay - 
ing a very low density (10 7_8 cm -3 , INicastro et al.|[l999l ). 
Alternatively, the outflowing gas may be "transient" , in the 
form of clouds or filaments passing through the line of sight. 
A variety of distinct ionisation stages, seemingly unrelated 
to the ionising flux, could be then accounted for. This ex- 
planation is strengthened if, contrary to what is observed 
for IZwl, components with very different column densi- 
ties were o b served at different tim es (e.g. IKraemer et al.l 
l200ll . 120051 ; ICostantini et al.1 l2000h . The outflow velocity 
may no t change in time, if the gas is also movi ng trans- 
versely (|Crenshaw et al.| [2003; Kr aemer et al JI2005D . Finally, 
the anti-correlation seen in Fig. [7] may be the result of a 
shielding of the gas which is thus experiencing a filtered ion- 
ising luminosity. Qualitatively, the (self-)shielding of the gas 
may happen away from our line of sight, and is not directly 
observabl e in a comple x gas environment as dep icted for ex- 
ample in IProeal (|2005l ): IProea Kallmanl l|2004h for the UV 
and X-ray emitting regions. 

4.3 The connection with the UV absorber 

Some of the gas components detected in the X-rays might 
also absorb in the UV band. In the H ST spectrum of I Zw 1, 
obtained in 1994 |Laor et al.lll997l ). the presence of Nv 
and C IV absorption features with an outflow velocity of 
— 1870 kms -1 , clearly highlighted the presence of a warm 
absorber. Unfortunately the column densities of the absorp- 
tion troughs were poorly determined in the HST spectrum: 
logAciv = 14.07 - 18.14 and log AW = 14.36 - 18.51. The 
derived total hydrogen column density was logA^H = 17 — 
20.77. The UV measurements (displayed with the dashed 
vertical lines) are compared in Fig. [8] with the values ob- 
tained in the X rays at different epochs. The ionic column 
densities of Nv and Civ, which are visible only in the UV 
band because of their low ionisation, were derived from the 
best-fit models of the X-ray absorbers. The a X-ray compo- 
nents have values of Anv and Nciv which are compatible 
with the UV column densities. Also the ASCA column den- 
sity measurements are consistent with the UV values, while 
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Figure 9. XMM-Newton velocity spectra for the 02a (up- 
per panel) and 05a (lower panel) components. The bin size is 
200kms —1 . The transitions chosen to produce the stack of the 
(unblended) line profiles for each absorber are also labelled. The 
vertical dashed line marks the blue-shift as measured by HST- 
FOS. Note that the redshift used in the HST analysis was 0.0608 
which implies a shift in the outflow velocity of only 90kms _1 , 
which is negligible compared to the errors associated to the X- 
ray measurements. 



the b components clearly do not produce enough of these 
low-ionisation ions. Interestingly, the UV flux in 1994 was 
about four times higher than in 02 and 05. This at least 
suggests that a low ionisation outflowing component may 
survive years, despite large changes in the UV flux (and 
presumably of the overall level of the SED). The UV to- 
tal hydrogen column density is slightly lower than for the 
X-ray absorbers (Fig. [SJ). The outflow velocity measured 
from the X-ray spectra is consistent with the value found 
in the HST-FOS spectrum. This is shown in Fig. [9] where 
the velocity spect ra for the 02a an d 05a components are 
compared (see e.g. iKaspi et al. 1 120021 , for details on the pro- 
cedure). The vertical dashed line marks the b lue-shift as 
meas ured from the UV Nv line (-1870 km s -1 : FLaor et all 
Il997l ). The outflow velocity of the UV-X-ray gas did not 
change in time, within uncertainties. This suggests that we 
are observing a gas not dramatically influenced by changes 
of the central source (e.g. changes in the acceleration of the 
flow, or strong variability of the source). Therefore the ob- 
served gas may be in the terminal-velocity portion of the 
flow, at a considerable distance from the central source, for 
example at a parsec scale, produced eit her by evaporation 
in a molecular torus ijKrolik Krissll200ll ). or by instabilities 
of the outer rim of the accretion disc, which in principle can 
extend to large distances. As a note of caution, the compari- 
son between the X-ray and the UV absorbers discussed here 
is b ased on non-simultane ous data, which may be mislead- 
ing |Crenshaw et al. ! 2003). In addition, the UV information 
given by HST gives only partial information on the SED, 
which would be important for an unambiguous comparison 
between the UV and X-ray bands. 



5 CONCLUSIONS 

We have presented the XMM-Newton low-energy X-ray data 
of I Zw 1, collected in 2005. In the absorbed spectrum we find 
evidence of a neutral phase, stable at least over the last ten 
years and probably associated with the host galaxy or with 
the interaction stream to the companion galaxy of I Zw 1 
(see also Paper I). In addition, two ionised absorption com- 
ponents with similar column densities of ~ few x 10 21 cm -2 
and ionisation parameters log£ ~ and 2.5 are detected. 
Comparing these results with a previous XMM-Newton ob- 
servation and archival ASCA and HST data we conclude: 

• The two ionised gas components detected in both 05 
and 02 cannot be in pressure equilibrium with each other. 
For the two components to coexist, a different mechanism, 
like magnetic confinement and/or gas shielding has to play 
a role. 

• The X-ray ionisation parameter £x is clearly anti- 
correlated with the X-ray luminosity. A viable explanation 
for this considers a gas organised in filaments crossing the 
line of sight. However the similarity of the column densi- 
ties of the a, b components and also of the gas detected 
in the ASCA spectrum, make this possibility not as clea r 
as for other cases (e.g. NGC 4151, iKraemer et al.l 120051 ). 
Other possibilities, like a non-equilibrium state of the ab- 
sorber, would need continuous monitoring to be rigorously 
tested. Finally, we may think that the observed gas might 
be shielded from the continuum source in the context of a 
complex and turbulent gas environment and experience a 
filtered ionising flux. 

• The gas producing the X-ray lower ionisation absorbers 
detected in 02 and 05 can also produce the lo w-ionisation 
speci es (N v and C iv) seen in the HST spectrum (|Laor et al.l 
1 1997T ) . as the column densities and outflow velocities are 
consistent. Interestingly however, the UV absorber was ob- 
served when the UV flux was about 4 times higher than 
the ones recorded by XMM-Newton-OM. This shows that a 
low ionisation component can survive for at least ten years, 
despite significant variations in the UV flux. The constant 
outflow velocity of the UV and the X-ray absorbers in time 
suggests also that we are observing the terminal velocity 
portion of the gas, probably located at a large distance from 
the source. The caveat is that the X-ray-UV comparison is 
performed with non- simultaneous observations, which might 
be misleading. 
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